A perturbational study of the lifetime of a Holstein polaron in the presence of weak 

disorder 
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Using the momentum average (MA) approximation, we find an analytical expression for the 
disorder-averaged Green's function of a Holstein polaron in a three-dimensional simple cubic lattice 
with random on-site energies. The on-site disorder is assumed to be weak compared to the kinetic 
energy of the polaron, and is treated perturbationally. Within this scheme, the states at the bottom 
of the polaron band are found to have an infinite lifetime, signaling a failure of perturbation theory 
at these energies. The higher-energy polaron states have a finite lifetime. We study this lifetime 
and the disorder-induced energy shift of these eigenstates for various strengths of disorder and 
electron-phonon coupling. We compare our findings to the predictions of Fermi's golden rule and 
the average T-matrix method, and find a significant quantitative discrepancy at strong electron- 
phonon coupling, where the polaron lifetime is much shorter than Fermi's golden rule prediction. 
We attribute this to the renormalization of the on-site potential by the electron-phonon coupling. 

PACS numbers: 71.38.-k,71.23.An,72.10.Di 



I. INTRODUCTION 

Studying the behavior of solid state systems under the 
simultaneous action of disorder and interactions is a sig- 
nificant challenge in condensed matter physics. Strong 
correlations in interacting systems often give rise to sharp 
quasiparticles. Scattering of such quasiparticlcs from 
weak disorder should just limit their lifetime. For strong 
disorder and no interactions, it is well understood^ that 
constructive interference of the backscattered waves can 
localize single particles such that they lose their itiner- 
ancy. If interactions are turned on, there is no consensus 
about the effect of disorder on the quasiparticles of inter- 
acting systems. 

For example, consider polarons, which are the focus 
of this work, and which are quasiparticles comprising a 
charge carrier plus a cloud of phonons describing the lat- 
tice distortions in the vicinity of the charge. If one views 
the polaron as a particle with a renormalized mass, large 
disorder should result in Anderson localization. How- 
ever, phonon-assisted hopping of carriers between local- 
ized states is well established as a conduction process in 
lightly doped semiconductors £ This suggests that, for 
suitable electron-phonon couplings, polarons may still be 
itinerant in a disorder potential that would localize par- 
ticles with the same effective mass. A possible explana- 
tion for this was offered recently in Refs. [H and [3, where 
the momentum average (MA) approximation was used to 
show that the electron-phonon coupling rcnormalizes the 
disorder potential in a strongly energy-dependent man- 
ner, so that the effective disorder seen by polarons can 
be drastically different from the bare disorder. 

References i and! focused on the effect of a single 
impurity potential on the polaron, although, as explained 
there, MA generalizes straightforwardly to other types of 
on-site potentials, including disordered ones. As such, it 
could be used to study numerically the effect of disorder, 



by generating results for various disorder realizations and 
analyzing their statistics. Such results have already been 
obtained for the Holstein^ polaron using the statistical 
dynamic mean- field theory (sDMFT)^l Given the quite 
different underlying approximations, it might be useful 
to check whether MA leads to similar results. 

A different approach, valid for weak disorder, is to 
use perturbation theory and perform the disorder av- 
erage analytically, similarly to the Born approximation 
widely employed for charge carriers in the absence of 
electron-phonon coupling^ Because localization cannot 
be described within a perturbational calculation, the po- 
laron eigenstates remain extended and self-averaging over 
all disorder realizations is appropriate. 

Here we follow the latter approach and calculate, using 
MA and for weak disorder, the disorder-averaged Green's 
function of the Holstein polaron and the resulting po- 
laron lifetime and energy shift. MA is an accurate an- 
alytical method originally developed for calculating the 
Green's function of Holstein polaron in clean systems^ 
and later extended to other types of coupling ! 10 : 11 MA 
is nonperturbative in the electron-phonon coupling as it 
sums all diagrams in the self-energy expansion, up to ex- 
ponentially small terms that are neglected. It also has a 
variational interpretation, in terms of the allowed struc- 
ture of the polaron cloudi 12 ' 13 MA can also be system- 
atically improved by increasing this variational spaced 
giving rise to the MA( ', MA' 1 ), etc., flavors which be- 
come more accurate but at an increased computational 
cost. Because here we focus only on the lowest-energy 
polaron states, which are already accurately described at 
the MA^ ' level, in the following we restrict ourselves to 
this flavor and call it MA for simplicity. 

The paper is organized as follows. In Sec. II, we 
present the generalization of MA to include disorder per- 
turbationally. Sec. Ill contains the results and their 
analysis, and Sect. IV has our conclusions. Various com- 
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putational details arc organized and presented in several 
appendixes. 



II. THE MODEL AND ITS SOLUTION 

The Hamiltonian for a single Holstein polaron in a lat- 
tice with random on-site energies is 



n = n d + v i-p h = n + v d + y c i_ ph , 



(i) 



where the non- interacting part of the Hamiltonian, H^, 
is divided into Vd = Cjcjci, describing the on-site dis- 
order potential experienced by the charge carrier, and 



H = -tj2(4cj+ H -c-)+nJ2 b i b * 



describing the kinetic energy of charge carrier plus the 
(optical) phonon energies (h = 1). The interaction part 



cl— ph 



E' 



(3) 



describes the Holstein coupling between the charge car- 
rier and phonons. As usual, q and bi arc annihilation op- 
erators for carrier and phonon, respectively, on the simple 
cubic lattice of lattice constant a, whose sites are indexed 
by i. The spin of the carrier is a trivial degree of freedom 
in this model, and we ignore it. Throughout this work 
we limit ourselves to the single polaron limit, i.e. there 
is a single charge carrier in the system. 

The on-site energies, {ei}, are taken from an uncorre- 
cted symmetric random distribution 



V({e i }) = U i V(e i ). 



(4) 



For Anderson- type disorder V(ti) has the customary 
form: 




if - A < ei < A 
otherwise, 



(5) 



while V(ti) = x8{ei — ea) + (1 — x)8(ei — ee) for a binary 
alloy, with x being the concentration of A-type atoms 
and energies shifted so that xtA + (1 — x)es — 0. 

Our aim is to calculate the Green's function of this Hol- 
stein polaron and average it analytically over all disorder 
configurations given by Eq. ([4]). From now on, we use 
an overbar to denote disorder-averaged quantities. The 



strength of disorder, a = 



is taken to be weak com- 



pared to polaron bandwidth in the clean system, so that 
it can be treated perturbationally. As a result, the pola- 
ronic picture remains valid but its lifetime is expected to 
become finite due to scattering from the disorder poten- 
tial Vd- It is precisely this disorder-induced lifetime that 
interests us. 

Since Vd is weak compared to other terms in the Hamil- 
tonian, we treat it as a perturbation. Dividing the Hamil- 
tonian as % = Hh + Vd , where "Hh is the Hamiltonian of 



the Holstein polaron in the clean lattice, we use Dyson's 
identity, G(u) = Gh(w) + G(w)VdGn(w), to relate the 
resolvent G(u>) of the system with disorder, to Gh(w) of 
the clean system. To the second order in Vd, we find 

G(u) w GhH + G H HVdG H H 
+G H (w)KiGH(w)VdG H (w). 

Because disorder breaks translational invariance, the 
eigenstates for any individual disorder realization are not 
labeled by the momentum k. However, averaging over 
all disorder configurations restores the translational in- 
variance and makes momentum a good quantum number 

<5k,k'G(k, w) and 



again. As a result, (0|ckG(cj)cj <; , |0) 



(2) we only need to calculate the diagonal matrix element: 



G(k, u) = G H (k, u) + ^(0|c k G H H4 Cl G H H4|0) 

i 

+ Y^(0\c k G K (^)4c i Gn(uj)c]c j Gn(u})ci\0}. 



hi 



Here, Gn(k, uS) = (0|ckG(cj)cj c |0) is the polaron Green's 
function in the clean system. For completeness, its MA 
solution is briefly reviewed in Appendix A. 

Since ei = for symmetric disorder, the first-order 
contribution vanishes (a finite average can be removed 
trivially by an overall shift of the energy). Because of 



uncorrclated disorder, 



tiki 



a Si j, and the 



disorder-averaged Green's function becomes 

G(k,w) = G H (k,w) 

+ g 2 ^(0|c k G H (u;)cJc i GHHctc i GH(w)4|0). 



The challenge is to use MA to calculate the matrix 
elements appearing in the second term to the same level 
of accuracy as Gn(k, w)J^ 

These matrix elements can be broken into products 
of generalized Green's functions by inserting identity op- 
erators, X, between the creation and annihilation oper- 
ators. Since the MA flavor we use here is equivalent 
with assuming that the phonon cloud only extends over 
one site^ at this level of accuracy it suffices to trun- 
cate! ps J^i n (l/ n} -)bj n \0)(0\b?, i.e., to ignore states with 
phonons at two or more sites (such states can be added 
systematically in higher flavors of MA). This leads to: 

G(k,w) =G H (k, W )+a 2 (0|c k G H Hc^"|0) 



i.i.s.n.m 



x (0|6fciG H (w)ct6t™|o)(0|^ Ci GH(w)c k |0)/(n!m!). 

This expression involves two sets of general- 
ized propagators, namely (0|ckGH(w)c^&}™|0) and 
(0|6J I c i G H (w)cj6j m |0). We now evaluate them. 
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First, as detailed in Appendix B, to the level of ac- 
curacy of MA, for n > 1 the first propagator van- 
ishes unless i = I. Therefore, we have to evaluate 
f£\ui) = (0|c k G H (w)c|6f l |0) for n > 1. Note that 
F^\u) is already known: F$(u) = (0|c k G H (^)4|0) = 
Gn(k, lo) exp(— ik • Rj)/\/iV, since momentum is a good 
quantum number in the clean system. Here, N — > oo is 
the number of sites in the system. The details of the cal- 
culation for n > 1, which is related to that of Gn(k, lo), 
are presented in Appendix B. The final result is: 



(6) 



where T n (uj) are easy to calculate products of continued 
fractions, see Eq. (|B1[) . 

Next, we calculate W nm (uj) = (0|^ l Cl G H (w)4^ m |0) 
to the same level of accuracy. Note that because of the 
invariance to translations in the clean system, this quan- 
tity is independent of i. The detailed derivation of these 
functions is presented in Appendix C. 

With these expressions in hand, the disorder averaged 
Green's function is, to second order in a: 



G(k,w) «G H (k,w)+ ( r 2 [G H (k,w)] 



n,m— 



r n (u)W n - m (u)T m (u) 



(7) 



To the same order, this identifies the sum in the previous 
equation as the disorder self-energy, 



use either A or a to characterize the disorder, as conve- 
nient, but we emphasize that any other type of disorder 
that has the same a would lead to the same answer within 
this perturbational approximation. To characterize the 
electron-phonon coupling strength, it is convenient to use 
the effective coupling A = g 2 /(6t£l). 

Once the Green's function is known, the energy broad- 
ening of a polaron state of momentum k, which is in- 
versely proportional to its lifetime, is the width of the 
low-energy peak in the spectral function, A(k, lo) = 
— i±mG(k, lo). This broadening measures the rate at 
which the polaron leaves that momentum state due to 
scattering from the impurity potential Vd- In a clean sys- 
tem, the polaron states are infinitely long lived, therefore 
the low-energy spectral weight is a Dirac delta function 
(in fact, a Lorentzian of width r) — >■ 0). Mathematically, 
this is a consequence of the fact that (in the absence of 
disorder) the polaron self-energy Ema(w) has a vanishing 
imaginary part for all energies inside the polaron band. 

Disorder-induced finite lifetime broadens the delta 
functions into Lorentzians. As a reference, we review first 
the case without electron-phonon coupling, A = 0. The 
only nonzero term in Eq. ([5]) corresponds to m = n = 0, 
therefore 



£ dis H = a 2 W°-°(uj) = a 2 g {uj), 



where 



flo 



V > Af 1,1 — <r. 



N t-r* lo - £k + if] 



(10) 



(11) 



S di »=a ^ — . (8) 

n,m— 



This is our main result. From a computational point 
of view, because the factorials in the denominator grow 
rapidly with increasing index, the infinite sums can be 
safely truncated at finite values for n and m. Cutoffs of 
20 proved sufficient for all cases we examined. 

Using Gn(k, lo) = 1/(lo — £k — £ma(w) + ii]), sec Ap- 
pendix A, we can finally write 



G(k,w) 



1 



lo - e k - StotM + iif 



(9) 



where the total self-energy is S to t(^) = Sma(w) + 
5jdis(w). This implicit summation gives a more accu- 
rate expression for the disorder-averaged Green's func- 
tion than Eq. ([7]), with which it agrees to 0(a 4 ). 



III. RESULTS 

We are now prepared to study the effect of weak disor- 
der on the polaron lifetime and energy shift. At this level 
of perturbation theory, disorder only enters through its 
standard deviation a. In the following, we assume Ander- 
son disorder of width 2A, for which a = A/-\/3. We will 



is the momentum-averaged free propagator. 

The resulting spectral weight has a peak of width T k 
centred at energy E^, found from the pole condition 



lo - £ k - EdisM + if] = 0, 



(12) 



where 



Because 7Y. 



r k 1 

T 2 



" and rj 



0^ 



- a is small for weak disorder, we ap- 
proximate Edis(w) w Edis(£ ; k) + C(cr 4 ). Using this in Eq. 
(|T2l gives _E k and t^ 1 as follows 

£k = £ k + RcSd is (^ k ) 
T k 1 = -Im£ dis (I? k ). 



(13) 



The first expression determines the energy shift com- 
pared to the electron energy in the clean system, £k- 
Since Hego(uj) is negative for lo < and positive for 
lo > 0, this implies a widening of the energy band in 
the presence of disorder. 

Using Eq. (fTU|) . the inverse lifetime becomes r^ 1 = 
— a 2 Imgo(E\ s ). However, Imgo(-E'k) is proportional to the 
total density of states (DOS) for the clean system, 



Imgo(-Ek) 



A 



') = -npo(Ek), 



so that: 



7n7 2 p (£k) = Tra 2 p (e k ) + 0{a 4 ). 



(14) 
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FIG. 1. (Color online) (a) Inverse polaron lifetime 1/rk vs 
its peak energy Ek, and (b) average DOS p(i^) vs uj, for a 
weak electron-phonon coupling and two values of the disorder 
A. The solid and dashed lines are the corresponding Fermi 
golden rule (FGR) and ATA results, respectively (see text for 
more details). Other parameters are fl = t,r)/t — 1CP 2 in (a) 
and r]/t = 5x 10" 3 in (b). 



FIG. 2. (Color online) Inverse lifetime of the polaron of mo- 
mentum k = (7r/8,0, 0) and A = 0.5 vs the strength of dis- 
order, A/t (full squares). The inverse lifetime for a free elec- 
tron with the same momentum is shown by triangles. Empty 
squares show l/(t*rk) vs A/t* for the polaron (for this A, 
t* — 0.88H). Other parameters are Q — t,r]/t = 10 -2 . 



The last equality is simply Fermi's golden rule. Since the 
density of states vanishes outside the bandwidth of the 
clean system, this result predicts infinite lifetime for all 
states with |£k| > 6i, and finite lifetime for all states in 
between. We will return to this point below. 

When the electron-phonon coupling is turned on the 
analysis is performed similarly, but now using the appro- 
priate total self-energy. The results are discussed next. 

We first consider weak electron-phonon coupling, A = 
0.5. In Fig. HJa) we plot the polaron inverse lifetime 
for states in the polaron band, for two different values 
of the disorder strength, A = 0.2t and OAt (squares and 
circles, respectively). These values are extracted from 
Lorentzian fits of the lowest peak in A(/k,co), using Eqs. 
(|13|) . The broadening rj was decreased until E^ and Tk 
converged to values independent of it. 

For this small A, the MA ground-state energy of the po- 
laron in the clean system is Ep^s = — 6.534i. The weak 
disorder does not shift the eigenstates significantly. In 
fact, as shown in Fig. QJb), the average density of states 
in the disordered system p{u) = — —lm ^ k G(k, uS) is 
nearly identical to that of the clean system, although the 
band becomes slightly broader with increasing A. The 
inverse lifetime vanishes below the clean system band- 
edge, Ep t cs, and above it increases like -JEu — Epjjs, 
which is the expected clean system DOS at the bottom 
of the band. This is very similar to the A = results, ex- 
cept for the renormalization of the DOS by the electron- 
phonon interactions. Indeed, if we think of the polaron 
as a simple quasiparticle whose density of states is p(ui) 
[renormalized from Po{lu) for a free electron], the inverse 
lifetimes we find at the bottom of the polaron band are in 
good agreement with those predicted by Fermi's golden 
rule (FGR), i.e., with -K<j 2 p{E\ s ) (see full lines). 



While the agreement between the two is good near the 
bottom of the band, it becomes systematically worse at 
higher energies. To verify that this amount of disorder 
is still sufficiently small so that the disagreement is not 
due to using perturbational results outside their valid- 
ity range, we also show average T-matrix (ATA) results 
(dashed lines). ATA is a simple way to treat disorder 
beyond the lowest order in perturbation theory, for a sys- 
tem with A = 0. We briefly discuss it in Appendix D, as 
well as how we extended it to finite A. ATA converges to 
ir(T 2 p(uj) in the limit of small c, therefore the agreement 
between FGR and ATA confirms that the contribution 
of higher order terms in a is indeed negligible. The dis- 
agreement with MA at higher energies is, therefore, not 
an artifact of using perturbation theory. 

The meaning of this disagreement at higher energies 
should, however, be treated with some caution. It is well 
knowrJ^ that this flavor of MA fails to reproduce the cor- 
rect polaron+one phonon continuum, which should start 
at Ep os + ^ (this problem is fixed by MA^ and higher 
flavors). One consequence is that MA overestimates the 
bandwidth of the polaron at weak couplings. Indeed, in 
Fig. [lja) we see that the polaron band extends well past 
EpfiS + ^- In other words, we know that at these higher 
energies MA is not accurate enough, so the results shown 
in Fig. []] should only be trusted close to the bottom of 
the band, where the agreement is good. 

The monotonic decrease of the polaron's inverse life- 
time with increasing disorder strength is shown in Fig. [2j 
for a polaron with momentum k = (7r/8, 0, 0) and A = 0.5 
(full squares). For comparison, also shown is the corre- 
sponding lifetime of a bare electron (A = 0, triangles) 
with the same momentum. Both curves show the ex- 
pected oc u increase predicted by Fermi's golden rule, 
but the polaron lifetime is somewhat shorter. The most 
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FIG. 3. (Color online). Top panels: l/(trk) vs E^/t for three 
levels of disorder: A/t = 0.05,0.1 and 0.15. The symbols 
shows the MA result for a strong coupling A = 1.2, while the 
full and dashed lines show Fermi's golden rule and the ATA 
pred ictions, respectively. Bottom panels: The average DOS 
p{uj) for that A (full line) and the DOS in the clean system, 
p(oj) (dashed line) vs u>. Parameters are fi = t,rj = 10~ 3 t. 



likely reason for this is the renormalization of the polaron 
mass by interactions. Indeed, if instead we plot l/(t*Tk) 
vs A/t* for the polaron (empty squares), the results are 
much closer to those of the free electron, especially for 
small values of the disorder. 

The conclusion, thus far, is that Fermi's golden rule 
agrees well with our results at energies where this fla- 
vor of MA can be trusted. In other words, at weak 
electron-phonon coupling, the effect of disorder can be 
quantitatively understood if we think of the polaron as a 
simple particle with a renormalized mass (or DOS), and 
use Fermi's golden rule. 

We now check whether this also holds true at strong 
electron-phonon coupling, for A = 1.2, where a robust 
small polaron appears in the clean system. The top pan- 
els in Fig. |3] show the polaron inverse lifetime vs. its 
energy for three levels of disorder (symbols), as well as 
the Fermi golden rule (full lines) and the ATA (dashed 
lines) predictions. The latter two are indistinguishable, 
confirming that these levels of disorder are indeed per- 
turbationally small. The bottom panels show the average 
density of states in the corresponding disordered systems, 
p{ijj) (full lines). For comparison, the polaron DOS in the 
absence of disorder, p(uj), is also shown (dashed line). 

Let us consider the DOS, first. As in the other cases, 
we see that with increasing disorder, the bandedgc shifts 
down to lower energies. However, the effect is quanti- 
tatively much more significant here than at weaker cou- 
plings because the polaron band is much narrower. This 
is seen in Fig. [4] where we show the same densities of 
states but over the full polaron band. One surprise is that 
the entire polaron band moves to lower energies with in- 
creasing disorder. This is different from what happens for 
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FIG. 4. (Color online) The same average DOS vs lj displayed 
in the lower panels of Fig. [3j but now shown for the entire 
polaron band. 



a simple particle, where the band broadens symmetrically 
on both sides. The different behavior at the upper edge 
is likely due to the difference in their spectra. While for 
a simple particle its band is the only feature in its spec- 
trum, the spectrum of the polaron is quite complicated, 
with many other features, such as a band associated with 
the second bound state, the polaron+one-phonon contin- 
uum, etc., lying above the polaron band.— With increas- 
ing disorder all these features should move toward lower 
energies. Level repulsion from higher-energy states would 
explain why the upper edge of the polaron band moves 
to lower energies. 

For the inverse lifetime we see that, as in the other 
cases, it vanishes for states with energy below the band- 
edge of the clean system, E'k < Ep^as ~ — 7.73£. Because 
the shift of the disorder-averaged DOS is now significant, 
this means that, for a quite large energy range at the 
bottom of the band, the polaron has an infinite lifetime 
despite the presence of disorder. We emphasize that this 
is qualitatively similar to the result for a simple parti- 
cle at the bottom of its band; the effect is simply quan- 
titatively more pronounced here. The meaning of this 
(un-physical) infinite lifetime for these low-energy states 
is discussed in the conclusions; briefly, we believe that 
it signals a failure of the perturbation theory at these 
energies. These low-energy states are most susceptible 
to localization, so the perturbational calculation and its 
predictions are suspect here. 

For higher-energy polaron states with E^ > Ep^cs, the 
lifetime in the presence of disorder becomes finite, as ex- 
pected. However, here the MA results disagrees quantita- 
tively with the FGR and ATA results at all energies. The 
latter two are nearly indistinguishable, suggesting again 
that these levels of disorder are small enough that pertur- 
bation theory should be valid. The disagreement cannot 
be blamed on MA in this case; at such strong couplings 
and correspondingly low energies, MA is extremely ac- 
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FIG. 5. (i) FGR approximation for the disorder-averaged 
Green's function of a carrier (double thin line), in terms of 
that of the free carrier (thin line) and uncorrelated disorder 
(dashed line) in the absence of electron-phonon coupling; (ii) 
FGR approximation for the disorder- averaged Green's func- 
tion of a polaron (double thick line), in terms of that of the 
clean polaron (thick line) ; (iii) clean polaron Green's function 
in terms of free carrier (thin lines) and phonon (curly lines) 
propagators; (iv) the first few terms in the MA approximation 
for the disorder averaged Green's function of a polaron. 



curate for the entire polaron bandi&H The disagreement 
is, therefore, meaningful. 

Its origin can be quite easily traced. If we explicitly 
separate the n = m = contribution in the disorder 
self-energy, Eq. (|5]) becomes 



SdisM = cr 2 g (uj - UmaM) 

n+m>0 



r n (Lo)W n > m (uj)T m (Lu) 



r I m I 



If the contribution of the terms with n + m > can 
be ignored, this result leads to Fermi's golden rule, since 
p(u) = — ^lmgo(uj — Sma(^))- The disagreement be- 
tween MA and FGR, then, comes from the contribution 
of the terms with n + m > 0. These terms cannot be ig- 
nored at strong electron-phonon coupling. Consider, for 

instance, F^\u) = (0|q k <V t ") c N n l >. which is pro- 
portional to T n (ui) and therefore is responsible for its 
appearance in £dis(w). If we Fourier transform to real 

times, iv? (t) is proportional to the amplitude of prob- 
ability that if an electron is injected in the system at 
some moment, at a time t later we find the electron in 
the presence of n bosons, all at the same site. At strong 
couplings, the electron dresses itself with a large phonon 
cloud to become a polaron, so the probability of finding 
it with many phonons in its vicinity should be consider- 



able, while the probability of finding the electron with- 
out any phonons (n = 0) is exponentially small. In the 
large A limit, the terms which are expected to contribute 
most are those with n rj g 2 /f2 2 , i.e., values close to the 
average number of phonons in the polaron cloud. In con- 
trast, for small A the phonon cloud is very fragile and, 
in fact, most of the time the electron is alone (resulting 
in a large quasiparticlc weight). This is why, for weak 
coupling, keeping only the n = m = term in the sum 
provides a good approximation. 

These considerations are illustrated diagramatically in 
Fig. [SJ Panel (i) shows the Born approximation for the 
disorder averaged Green's function of a carrier, which 
leads to Fermi golden's rule expression for the lifetime, 
as already discussed. Panel (ii) shows its equivalent for 
the disorder averaged Green's function of the polaron; 
as discussed above, this is equivalent with keeping only 
the n = m = term in the disorder self-energy, Eq. @. 
Since each clean polaron propagator starts and ends with 
a free carrier propagator, this approximation means that 
the electron can scatter on disorder only in the absence 
of phonons; this is why this approximation fails at large 
electron-phonon coupling, where a large phonon cloud 
forms. In contrast, the full MA expression includes di- 
agrams such as shown in panel (iv), where phonon and 
disorder lines cross. One can think of these as leading 
to an effective renormalization of the disorder strength, 
especially since these diagrams are very similar to those 
which result in the renormalization of a single impurity 
potential, discussed in Refs. and [3. 

Our results in Fig. [3] reveal that this renormalization 
of the disorder (considerably) lowers the polaron lifetime. 
Thus, the strong electron-phonon coupling effectively en- 
hances the strength of the disorder potential. Why this is 
so is not easy to infer from the results for a single impu- 
rity potential, where renormalization can either increase 
or decrease the bare potential, and even change its sign, 
and all these behaviors are seen at different energies (re- 
tardation effects are very significant). In any event, Fig. 
[3] shows that this renormalization is quantitatively sig- 
nificant. 

A more detailed picture of the evolution of Tk and £k 
with disorder A at strong coupling A is shown in Fig. [6] 
for two momenta ki = (2tt/9, 0,0) and k 2 = {n/6, 0,0), 
which have free electron energies £ki ~ — 5.5t and £k 2 ~ 
—5.7t, respectively. In Fig. [Hla), we trace their inverse 
lifetimes as a function of disorder. At small A both these 
states are well above Ec,ps, an d their scattering rates 
increase monotonically with A, as one would expect on 
general grounds. However, the inverse lifetimes reach a 
maximum after which they begin to decrease fast and 
eventually vanish. The value of A where they vanish is 
the disorder at which that eigenstate has £k = Ep } cs- 
For larger disorder, this eigenstate moves below the free 
polaron bandedge, and its lifetime becomes infinite. This 
is more clearly shown by Fig. [^b), where the inverse life- 
times are plotted vs. the corresponding cigenencrgy £k 
for these two momenta, as disorder is increased. This 
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l/tx, 



2x10 



1x10 




FIG. 6. (Color online) (a) Inverse lifetime vs disorder, and 
(b) inverse lifetime vs energy E^, as disorder is turned on, 
for two momenta ki = (2-k/9, 0, 0) and k2 = (7r/6, 0, 0), for a 
polaron with A = 1.2, Q, = t, rj = 10~ 3 t. The asterisk in panel 
(b) marks the clean polaron GS energy in the clean system, 
Ep,gs, for these parameters. See text for more details. 



confirms that the scattering rates for both polaron states 
vanish when their energy drops below Ep^qs, whose lo- 
cation is marked by the asterisk. The value of A where 
this happens depends on how far above Ep^cs was the 
energy E^ of this polaron, in the limit A —> 0. 

Figure [3] showed that using the FGR estimate, i.e., 
t^T 1 = 7rer 2 p(Ek), is quantitatively wrong. We can also 
compare the polaron lifetime, where finite, with that of 
a free particle of rcnormalized mass, similar to the com- 
parison in Fig. [2] This is shown in Fig. where we 
compare l/(t*Tk) vs A/t* for the polaron, with the in- 
verse lifetime of a free electron with the same momentum. 
While roughly quadratic dependence is observed for the 



l/t*x 



0.002 



0.001 



□-□polaroii, h=l. 2, t*=0.148t 
free electron, X=0, t*=t 




FIG. 7. (color online) Same data as shown in Fig. |6ja) but 
with rescaled axes; l/(t*Tk) vs A/t* for momentum k2 = 
(27r/9, 0, 0) (empty squares) is compared with the free electron 
lifetime (triangles) at the same momentum, for A = 1.2 where 
t* — 0.148£. Other parameters are Q — t and rj/t — 1CP' 5 . 



polaron at small disorder, the coefficient is quite different 
from that for the free electron. At higher disorder, the 
disagreement is even worse. 

This shows that for intermediate and large electron- 
phonon coupling, where a heavy small polaron forms, 
its lifetime in the presence of disorder is not described 
quantitatively by the predictions corresponding to a sim- 
ple particle with renormalizcd mass. The polaron has 
an internal structure which manifests itself in significant 
corrections to Fermi's golden rule even for weak disor- 
der. The scattering of the electron in the presence of its 
phonon cloud is quite different from that of a simple par- 
ticle of the same effective mass, but without a cloud 4^ 



IV. SUMMARY AND CONCLUSIONS 

Using MA to deal with the electron-phonon coupling 
and perturbation theory to deal with the weak disor- 
der, we derived an expression for the disorder-averaged 
Green's function of the Holstein polaron in a simple cu- 
bic lattice with random on-site energies. This allowed us 
to find an analytic expression for the lowest-order contri- 
bution from disorder to the polaron self-energy. 

The disorder-averaged spectral weight was used to ex- 
tract the lifetime and energy shift of various polaron 
states. For weak electron-phonon coupling, we found that 
the MA results are in reasonable quantitative agreement 
with those predicted by Fermi's golden rule for a free 
particle with an appropriately rcnormalized mass. 

At intermediate and larger electron-phonon coupling 
where a small polaron forms, however, the MA results 
quantitatively disagree with Fermi's golden rule estimate 
everywhere the lifetime is finite and for all levels of disor- 
der. The reason for this is the fact that the scattering of 
the electron in the presence of its (robust) phonon cloud 
is quite different from the scattering of a simple parti- 
cle with renormalized mass. This is the same physics 
that leads to a significant renormalization of the disor- 
der potential seen by a polaron as compared to the bare 
disorder .^ii This demonstrates that, in the small polaron 
limit, it is wrong to assume that the only effect of the 
polaron cloud is to renormalize the polaron's mass. 

It is important to note that this calculation is only 
valid for weak disorder. It is based on perturbation the- 
ory, and in principle it can be improved by going to 
higher orders along the same lines we used to calculate 
the lowest-order contribution. However, one should re- 
member that as disorder becomes stronger, Anderson lo- 
calization will eventually occur, and that this cannot be 
captured within perturbation theory Also, once disor- 
der is large enough to lead to localization, the disorder- 
averaged Green's function loses its meaning and useful- 
ness. Instead, here the signature of localization becomes 
manifest in the distribution of various quantities such as 
the local density of states, not in their average value. 

A surprise, at least at first sight, is the fact that this 
calculation predicts an infinite lifetime for a range of en- 
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crgies at the bottom of the polaron band. This interval 
can include a significant fraction of the polaron states, es- 
pecially at stronger electron-phonon coupling and larger 
disorder. As we already mentioned, this is in fact similar 
to what happens for a free particle, which also is pre- 
dicted, within this level of perturbation theory, to have 
an infinite lifetime for all momenta for which |i? k | > 6t. 
The difference is only quantitative: the energy shift for 
a free particle is tiny compared with its 12t bandwidth, 
whereas for a small polaron this shift can be comparable 
with its significantly narrower bandwidth even for rather 
weak disorder. 

A likely reason for this can be inferred from the fact 
that, for a free particle, states at the band edge become 
localized immediately upon introduction of disorder. In 
other words, we already know that there is a finite range 
of energies (which, for weak disorder, falls outside the 
free particle bandwidth) where treating disorder pertur- 
bationally and calculating the disorder-averaged Green's 
function is meaningless. It is then reasonable to conclude 
that the states for which this pcrturbational scheme pre- 
dicts infinite lifetimes are, in fact, already localized. If 
this is correct and generalizes to the polaron case, it sug- 
gests that, unlike for a free particle, for a polaron localiza- 
tion sets in differently at the lower vs the upper polaron 
bandcdge. We have already speculated that this differ- 
ence may be due to the influence of the higher-energy 
states that exist in the polaron spectrum. Confirmation 
of these conclusions will require a study going beyond a 
pcrturbational treatment of disorder. 
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Appendix A: MA solution for the clean system 

Using the same notation as in the main part, the Hol- 
stein Hamiltonian for a clean system is 



Hh = Ho + Vel-ph- 



(Al) 



To calculate G H (k,w) = (0|c k G H (u;)4|0), we use 
Dyson's identity Gh(w) = Go(w) + Gn(uj)V c i- p hGo(u}). 
The propagators for Ho are (0|c k &™Go(w)c k ,6j m |0) = 
n!5ri,m^i,i^k,k'Go(k,a; - nfi), where 



G (k,w) 

We find 

G H (k,w) = G (k,w) 



id - £k + ^ 



N 



(A2) 



where we introduced the generalized propagators 
= (0|c k G H Hc!&hO). Note that F$(u) = 
GH~(k, uj) exp(— ik ■ Rj)/\/iV. Using the Dyson identity 
again, we find that for n > 1 



F, 



+ gGo{i, i,uj — nO) 



g^G (j,i,uj-nn){0\c k GK(u})c]b]bt n \0) 



nF, 



(n-l) 



M + f£ +1) M 



Here 



G (j,i,w) 



k 



k(Rj-Ri 

N 



-G (k,w) 



is the free propagator in real space (the sum is over the 
Brillouin zone). This propagator decays exponentially 
with the distance |Rj — R,| for energies outside the free 
particle continuum, |cj| > 6i. Since we are interested in 
energies uj — nfl ~ Epqs — nil, where the polaron GS 
energy Ep qs < —6t, all these propagators become expo- 
nentially small for j ^ i. It is therefore a reasonable first 
approximation to ignore j ^ i terms in the EOM written 
above. This is what the MA approximation does (higher 
flavors include j^i terms in a certain progression^). 
Within MA, then, we have for any n > 1 



^IwM = -950 (^ - nil) 



nF 



(n-X) 



(A3) 

where go{ui) = Go(i,i,w); see Eq. (| 1 1[) . This recurrence 
relation is solved in terms of continued fractions^ 



(u)=A n {u>)Fi n r 1 \u) 



for any n > 1, where 
A n {u)) 



nggo{uJ - nfi) 



I - gg (u> - nCl)A n+ x(uj) 



Finally, using F£/(w) = Ai(w)F^V(w) in Eq 
to the MA solution: 



(A4) 



(A5) 



2h leads 



G H (k,w) 



where 



to — e k - SmaM + ir) 
Sma(w) = gAx{u). 
Appendix B: MA solution for F^™\uj) 



(A6) 



As discussed in the text, we need to calculate the 
propagators: (0 1 c k G H (w) c| b \ n |0). Within the MA ap- 
proximation (see above), this propagator is set to zero 
for all i ^ s, because it is proportional to Go{s,i,u> — 
nfl). The only finite value is for i = s, in which 
case (0|ckG H (w)ct&t re |0) = F$(w). These propaga- 
tors have already been calculated above, Fj£ \ui) = 
4»i& -1) M = • • • = T n (u)F^\u), where ' 

r n (w) = Ki^K-iiu) ■ ■ ■ Ai{u). (BI) 
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Appendix C: MA expression for W n,m (ui) 



Here we calculate the remaining needed propaga- 



tors, W nm {uj) = (0\b^CiGu(^)clb, 



Ut™l 



W mn {uf). 



Because of this symmetry, we only need to find 
W nm (ui) for to < n. Note that we already know 
W 00 (lo) = (0| Ci GhM4|0) = ^E k G H (k, W ) = g (u- 
EmaM), and also W n0 (oj) = (0|6?c,-6h(w)4|0) = 



F M ( w )l»7->- »7 = r n( w )9o(^ - SmaM)- 

For any m > 1, writing the EOM for W" m (o;) within 
the MA approximation (i.e., not allowing the electron to 
change its site), leads to 

W n ' m (ui) = m\g (uj — mQ)S nm + ggo(w — mCl) 

x{mW n > m - 1 (uj)+W n < m+1 (uj)}. (CI) 



For m > n, the delta function vanishes and Eq. (|C1|) is 
identical to Eq. (|A3|) . hence 



W n ' m (w) = A m {u)W r ' 



(C2) 



In particular, this gives W n ' n+1 {uj) = A n+1 (cj)W n > n (u) . 
Using this in Eq. (|C1|) with n = to relates W n ' n {w) to 



W n ' n (tu) = A n ((j)W n ' n ~ 1 (cj) + 



(n-l)l 
9 



A n (u). (C3) 



This is taken together with the EOM for 1 < m < n— 1 
W n ' m (u)) = ggo(u-mn) x [roW 1 '™ -1 ^) + W n ' m+1 (u)} 

to give a system of n equations with n unknowns 
W n ' m (u)), m = l,..,n [W n '°(u)) is known; see above]. 
This can be solved in many ways, including direct numer- 
ical solution. A nicer approach is to use the linearity of 
this system of equations to split it into two different sys- 
tems, one which has only ^-^-A n (u>) and one which has 
only ggo(uj — fl)W n '°(u>) as inhomogeneous parts. These 



W n ' m (co) = T m (u)W nfi {u) +T m (cu) 



(n-l)L4„H 
g[l - A n {u)B n (u)Y 



where T m (u) = B m+1 (Lu)B m+2 {u) ■ ■ ■ B n {oj) for to < n 
while r„(w) = 1, and 



B m +i{w) 



99a (w - mfi) 



1 — {tn — l)ggo(ui — m£l)B m (uj) 
are continued fractions ending at i?2(w) = ggoi^ ~ ^)- 

Appendix D: ATA self-energy 

AT A 15 i 16 relates the disorder part of the self-energy of 
a single particle to the disorder average of its transfer- 
matrix through a single impurity with on-site energy e: 



Sata(w) 



t 



(Dl) 



where t = e/(l — ego(uj)) is the sum over all single im- 
purity scattering contributions. go(w) is the momentum 
average of the free particle propagator, see Eq. (fTT|) . 
For Anderson-type disorder we find 



I 
2A 



ede 
1 - eg (oj) 
1 



hi- 



- Ago(u) 

'sb(w) T 2A. g 2( w ) 111 1 - A ff oM' 

Expanding to lowest order regains the perturbational 

result, Eq. E A ta(w) « ^flo(w) = ^g^)- Dif- 

ferences between ATA and FGR show that disorder is 
so large that multiple scattering processes off the same 
impurity cannot be ignored anymore. 

To extend ATA to the Holstcin model, we note that 
difference between the MA clean polaron's Green's func- 
tion and that of the free electron is the appearance of 
Ema(w). This simply modifies go(u>) — > go(tu — Ema(w)) 
in all the above equations. As discussed above, this ap- 
proximation implies that there is no crossing between 
phonon lines and scattering lines, in diagramatic terms. 
Our results show that this is a bad approximation at 
larger electron-phonon coupling. 
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